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We present an algorithm for improving subsequent parton shower emissions by full SU(3) 
colour correlations in the framework of a dipole-type shower. As a proof of concept, we 
present results from the first implementation of such an algorithm for a final state shower. 

1 Introduction and motivation 

Parton showers and event generators are indispensable tools for predicting and understanding 
collider results [TJ[2[3]. Considering their importance for interpreting LHC results, it is essential 
to have a good understanding of their approximations and limitations. These simulations have 
up to now all been based on QCD as an SU(N C ) gauge theory in the limit of large N c . For 
N c = 3, this approximation seems to work remarkably well despite the fact that 1/N C = 1/3, 
in the general case, and 1/N% = 1/9 in most cases are not truly small parameters. Otherwise, 
significant deviations from parton shower predictions as compared to experimentally measured 
observables would have already hinted towards a severe underestimate of colour suppressed 
terms. 

Including colour suppressed terms in parton shower simulations has so far been an unex- 
plored continent, and investigating effects caused by colour correlations beyond the large- N c 
limit is mandatory in the age of ever improving simulations, particularly when including higher- 
order QCD corrections. Especially when considering the matching of parton showers to NLO 
QCD corrections, subleading-iV c improved parton showers provide a valuable input in making 
these matchings more precise such that the matching conditions are indeed satisfied exactly, 
and not only modulo colour suppressed terms Q 

We here present an approach to subleading colour contributions [5] which is simple in the 
sense that it fits very well into the framework of existing Monte Carlo event generators. We 
note that this is not the end of the story, as for the general case an evolution at amplitude level 
would have to be considered. Our approach of colour matrix element corrections is a first step 
towards quantifying the size of the expected effects. 

2 From dipole factorization to dipole showering 

Dipole factorization, [BJ, states that the behaviour of QCD tree-level matrix elements squared 
in any singly unresolved limit involving two partons i,j (i.e. whenever i and j become collinear 

1 In this context an independent approach, considering only one emission, has been presented in [4]. 
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or one of them soft), can be cast into the form 

\M n+ x(...,Pi, ...,pj, ...,p k , ■•■)| 2 ~ 

— (M n {...,p {p ...,p k ,...)\V ihk (p l ,p 3 ,p k )\M n {...,p- J ,...,p k ,...)) , (1) 

where \A4 n ) - which is a vector in the space of helicity and colour configurations - denotes the 
amplitude for an n-parton final state. Here an emitter ij undergoes splitting to two partons 
i and j in the presence of a spectator k which absorbs the longitudinal recoil of the splitting, 
k — > k. This factorization formula, which is well established to provide a subtraction scheme for 
NLO calculations, can actually be used to derive a dipole-type shower algorithm [7]. Results of 
an implementation have been reported in [5], and similar approaches have been considered in 
[Hllini- In these cases, the colour correlations present in Vy * are approximated in the large- 
N c limit, while keeping the colour factor for gluon emission off quarks, Cf — T 2 ,. , exact. In 
turn, chains of colour connected dipoles evolve through subsequent emissions generating more 
dipoles in a chain or leading to a breakup of the chain in case a gluon splits into a qq pair until 
eventually the transverse momentum of potential emissions is below an infrared cutoff in the 
region of one GeV. 



3 Colour matrix element corrections 

The dipole factorization formula eq. [l] implies a factorization at the level of cross sections; here, 
the differential cross section for n + 1 partons factorizes into the cross section for producing n 
partons times a radiation density as the sum over all dipole configurations ij, k, which undergo 
radiation, ij, k — » i, j, k: 



, ^ n+1 (pl,z) -1 (M n \T f] -T k \M n ) 
Tfi \M n \ 



dp ijik (pi, z) = % fe ( P i, z) —^ ± ' >< :l (2) 

'ij 



Here, we have used the spin-averaged version of the dipole kernels including the product of 
colour charges encoding the colour correlations, V,,-fc = — fcTr- -T^/T 2 -., and dcj) k denotes 
the fc-parton phase space. In the large-A^ limit, this formula yields the basis for the dipole 
shower considered so far: — Tr. ■ T k /T 2 r. —> S(ij, k colour connected)/ (1 +6^), where 8q = 1(0) 

for ij = g{q/q). To obtain an algorithm which instead keeps the full colour correlations, we 
do not consider this approximation but keep the second factor in eq. [2] exactly. Owing to 
the similarity of matrix element corrections present in parton showers so far, we refer to this 
improvement as 'colour matrix element corrections'. 

Eq. [2] describes how a single emission incorporates colour correlations. Indeed, for the 
first emission off the hard subprocess, \M n ) is known, though it has to be recalculated after 
each subsequent emission to define the colour matrix element correction for the next emission. 
Instead of directly calculating the next amplitude, which would only be possible if we had 
derived splitting amplitudes in the singly unresolved limits, we observe that 

\M n \ 2 = M\S n M n = Tr (S n x M n M\) (3) 

and 

(M n \T rj ■ T- k \M n ) = Tr [S n+1 x T kn M n MiT} 3 n ) , (4) 
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where we have chosen a definite basis {|a)} for the colour space, \M n ) = 2 a =i c n,a\ct n ) <^ 
M-n = {c-n.ii •••) Cn,d n ) T an d introduced the scalar product matrix S n = {(a n \/3 n )} as well as 
matrix representations of the colour charge operators for n partons, Tj — > T i n . 

These representations then imply that we can work with an amplitude matrix M n as the 
fundamental object, 



where the initial matrix for the hard subprocess is given by M„ hard = .M nhard .Mjj h . 

4 Technicalities 

Having outlined the principle of the algorithm for including colour correlations for subsequent 
parton shower emissions, two major technical issues have to be addressed: On the one hand, 
an general treatment of the colour basis for an arbitrary number of partons is required. On 
the other hand, sampling from the probability (Sudakov-type) density driving the next parton 
shower emission has to be generalized to the case of non-positive splitting rates as typically 
encountered for l/N c suppressed contributions. 

For the first task, we have implemented a C++ library ColorFull 11 implementing the trace 
bases of colour space, [T2]. This library is interfaced to the Matchbox framework presented in 
[8]. The colour matrix element corrections calculated in this part of the simulations, are inserted 
as correction weights into an existing dipole shower implementation, which uses the ExSample 
library |13j to sample Sudakov-type densities derived from the absolute value of the colour- 
corrected splitting rates. For the second task, we then employ the interleaved competition/veto 
algorithm outlined in |14j to arrive at events distributed according to the desired density (note 
that the sum of all splitting rates approximates a squared matrix element and is thus positive). 

5 Results 

As a proof-of-concept, we present results from the subleading-iV c improved parton shower for 
final state radiation, more precisely considering e + e~ — > qq at LEP1 energies. We compare 
three different approximations: 'full' colour correlations, the 'shower' approximation where the 
TVTj/Tf are taken in the large--/V c limit, and a 'strict' large-A'c- approximation, where also 
Cf ~ CU/2. Interestingly, the 'shower' approximation does not exactly reproduce the shower 
implementation; from four partons onwards, this evolution is sensitive to the emission history, 
though matches the shower implementation if one sequence of dipoles has dominated. Indeed, 
the differences between the 'shower' approximation and the default shower implementation are 
at the per-mille level as one would expect from strong ordering in the emission history. 

In the results presented here, we include up to six improved shower emissions, g — > qq 
splittings are neglected, as there are no associated colour correlations; we also do not include 
hadronization. For event shapes and jet rates we find small subleading-iV c effects when consid- 
ering the shower approximation; for tailored observables, probing the dynamics of soft radiation 
with respect to a hard subsystem of the event, larger effects are seen. The strict approxima- 
tion shows larger deviations. This fact can mainly be attributed to the change in the Sudakov 
exponent for gluon emission off a quark. A few sample results are shown in figure [I] 
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Thrust , T = 1 - T 



average rapidity w.r.t. 773 




Figure 1: The thrust distribution (left) and the average rapidity w.r.t the thrust axis defined 
by the three hardest partons using the different approximations. 

6 Conclusions 

We have presented the first implementation of a subleading-./V c improved parton showers. For 
e + e~ — s- jets small effects are seen except for very special observables. The technical issues 
associated with the implementation, particularly the treatment of the colour basis and the 
presence of negative splitting kernels will serve as input for related and future work; we also 
anticipate that larger effects can be seen in hadron collisions, e.g. pp — > jets. 
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